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Using a combination of analytical theory, Monte Carlo simulations, and three dimensional self- 
consistent field calculations, we study the equilibrium properties and the switching behavior of 
adsorption-active polymer chains included in a homopolymer brush. The switching transition is 
driven by a conformational change of a small fraction of minority chains, which are attracted by the 
substrate. Depending on the strength of the attractive interaction, the minority chains assume one of 
two states: An exposed state characterized by a stem-crown-like conformation, and an adsorbed state 
characterized by a flat two-dimensional structure. Comparing the Monte Carlo simulations, which 
use an Edwards-type Hamiltonian with density dependent interactions, with the predictions from 
self-consistent-field theory based on the same Hamiltonian, we find that thermal density fluctuations 
affect the system in two different ways. First, they renormalize the excluded volume interaction 
parameter Vbare inside the brush. The properties of the brushes can be reproduced by self-consistent 
field theory if one replaces iWe by an effective parameter v e s, where the ratio of second virial 
coefficients _B e ff/-Bbare depends on the range of monomer interactions, but not on the grafting density, 
the chain length, and Wbare- Second, density fluctuations affect the conformations of chains at the 
brush surface and have a favorable effect on the characteristics of the switching transition: In the 
interesting regime where the transition is sharp, they reduce the free energy barrier between the two 
states significantly. The scaling behavior of various quantities is also analyzed and compared with 
analytical predictions. 


I. INTRODUCTION 

Polymer brushes are organic surface layers formed by 
polymer chains that are grafted at one end to a sub¬ 
strate 0 , 1 . Since the thickness of brush layers is typi¬ 
cally in the nanometer range, they are interesting for the 
design of functional surfaces [3H5J with applications in 
a wide variety of areas ranging from colloidal stabiliza¬ 
tion 1, lubrication [fi), controlled friction and adhesion 

& anti-fouling [tj, biocompatibility [l(|, drug delivery 
, and smart stimuli-responsive materials [4|. In this 
respect, multicomponent polymer brushes are particu¬ 
larly promising [4|, j5| . If the brush chains are covalently 
bound to the substrate, they cannot phase separate on a 
global scale, but they can still develop structure on the 
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nanoscale. The resulting brush morphologies are con¬ 
trolled both by the intrinsic properties of the brush, e.g., 
chemical properties of the chains (compatibility or in¬ 
compatibility), the grafting densities, the chain lengths, 
and by the environment-related parameters, such as sol¬ 
vent selectivity, substrate preference, temperature, and 
the pH. Due to the ability of polymer brushes to se¬ 
lectively respond to environmental stimuli, they can be 
used to design materials that can reversibly switch/tune 
their surface properties, e .g., with respect to wettabil¬ 
ity 01 p ermeability [l4j j, friction jl5|] . and optical 
properties jl6j |. 

Stimuli-induced phase separation provides the basic 
mechanism for the change in the brush surface proper¬ 
ties depending on which of the two microphases forms 
the outer part of the brush. A typical example for such 
a morphology-related switchable surface is a mixed poly¬ 
mer brush with equal amounts of hydrophobic and hy¬ 
drophilic polymers grafted on a substrate. By treatment 
of different solvents, the surface composition can change 
and then the wettability of this material switches [121. 113| . 
Morphology change necessarily involves slow highly coop- 
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erative chain dynamics and therefore a typical response 
time turns out to be in the range of several minutes or 
larger [jj| . 

In a recent letter, we have proposed a new class of 
brush-based switches 0 , which rely on a radical confor¬ 
mational change of individual adsorption-active minority 
chains in an otherwise inert brush. 

The selective adsorption driving the transition may 
arise, e.g., from electrostatic interactions or hydrogen 
bonding between active groups in the minority chains and 
the substrate. These active groups can serve as respon¬ 
sive sensors that trigger the switching transition. The 
most radical change is associated with the end group of 
the minority chain. In the adsorbed state it resides in 
close proximity of the substrate deeply buried within the 
brush, while in the other state it is exposed to the envi¬ 
ronment. Thus the transition could potentially promote 
a specific immune-like response. Chemically or biologi¬ 
cally active groups can be attached to the free end of the 
minority chain to serve as practically useful sensors. 

Based on theoretical arguments and simple one¬ 
dimensional mean-field calculations, we demonstrated 
that even a small chain length increment of about 10 
% produces a sharp transition from the exposed state to 
the adsorbed state. The transition time is expected to be 
very short since the free energy barrier between these two 
thermodynamically stable states is about several /c^T. A 
further increase in the minority chain length will lead to 
sharper transitions, but at the same time to longer transi¬ 
tion times due to higher barriers. The stron g re sponse of 
the chain conformation of a minority chain [181 ] to small 
variations of its length is related to the recently reported 
“surface instabilities” in polymer brushes [19|, 120 ] , which 
can be used to sense solvent quality 00. In our work, 
we proposed to use single chains as switches triggered 
by a change in substrate-polymer interaction. One ma¬ 
jor benefit of this switch is that it does not involve co¬ 
operative rearrangements of many chains, therefore the 
switching transition is fast as compared to the existing 
examples in mixed brushes. 

The arguments presented in Ref. 0 rely on several 
assumptions. First, the presence of the minority chains 
was taken to have no effect on the surrounding poly¬ 
mer brush. Second, thermal fluctuations of the majority 
brush component were disregarded. 

In the current paper, we present extensive Monte Carlo 
(MC) simulations of a coarse-grained model for polymer 
brushes with a single immersed adsorption-active minor¬ 
ity chain. To assess the influence of fluctuations sep¬ 
arately, we have also performed three-dimensional self- 
consistent field (SCF) calculations for comparison. We 
investigate in detail both lateral and longitudinal char¬ 
acteristics of chain conformations and the various char¬ 
acteristics of the transition. 

The remainder of the paper is organized as follows: 
Sec.[TT]outlines the MC model and the simulation scheme. 
Sec. M presents the MC results for a homogeneous 
monodisperse brush and comparison with the established 


analytical SCF theory, and analyzes the dependence of 
the renormalized virial coefficient on the effective inter¬ 
action range. Sec. IIVI presents a more detailed derivation 
of the theory for the adsorption-active minority chain 
sketched in Ref. 0. In Sec. 0 the main results of the 
MC simulations are presented and discussed. This sec¬ 
tion also compares MC and 3-d SCF calculations with 
the emphasis on fluctuation effects. A general discussion 
is given in the final section IVII 


II. MC SIMULATION MODEL 

In the simulation community, a variety of numeri¬ 
cal methods have been developed and used to inves¬ 
tigate polymer brush systems, ranging from molecular 
dynamics 0], dissipative particle dynamics [0, and 
Monte Carlo (MC) simulations [251 - 12911 to numerical self- 
consistent field (SCF) calculations 13(11351 . The present 
study mainly relies on MC simulations, while for com¬ 
parison we also present results obtained from 3-d SCF 
theory. The MC simulation model and scheme are de¬ 
scribed as follows, and the detailed description of SCF 
method is shown in the Appendix. 

In the MC simulations, we adopt a coarse-grained off- 
lattice model first proposed by Laradji et al [26[. In this 
approach, a particle-based representation of the poly¬ 
mers is combined with an Edwards type Hamiltonian 
[HI Hi} j which defines non-bonded interactions in terms 
of local monomer densities. Compared to the more com¬ 
monly used coarse-grained polymer models with hard¬ 
core monomer interactions, the Laradji model has two 
advantages: First, it can be simulated very efficiently, 
since the explicit evaluation of pair interactions is often 
the most time consuming part in a simulation. Second, 
it uses soft potentials, hence equilibration times are com¬ 
paratively short. The model does not account for packing 
effects and for topological constraints (which may restrict 
the conformational phase space for strongly adsorbed, 
quasi-two dimensional polymers). However, these are not 
in the focus of the present study. 

Specifically, the model system is a monodisperse brush 
in a good solvent containing a single minority chain in 
a volume V = L x ■ L y ■ L z . We use periodic boundaries 
along x and y directions, while impenetrable boundary 
walls are placed at 2 = 0 and z = L~. Polymer chains 
are modelled by the discrete Gaussian bead-spring model 
with spring constant 3 , where a is the statistical bond 
length, ks is the Boltzmann constant, and T the temper¬ 
ature. We will use a as the basic length unit, and as 
the energy unit. All chains are attached to a substrate 
placed at zq at one end, where zq is chosen smaller than 
a. The impenetrable wall at z = 0 exerts an attractive 
potential to all minority chain beads with strength £ in 
a range of 0 ^ z ^ a. There is no explicit solvent in the 
model and all the relevant interactions in the good sol¬ 
vent case are represented through an effective excluded 
volume potential between monomers (beads). The Ed- 
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wards type Hamiltonian of the system reads 

o n b Nb - 1 

^=2?EE (Raj — Raj'-r) 2 

a=l j=l 
JV-1 

+ ^E( R 0i- R 0d-i) 2 
i=i 

+ J dr p 2 ^ + J drU ads (r)p m (r), ( 1 ) 

where /3 = 1/ksT, R aj - denotes the position of the j- th 
bead in the a-th chain (the minority chain has the in¬ 
dex a = 0), nb is the total number of brush chains, p m 
denotes the local density of minority chain monomers, 
Pm = J2j S(r - R 0 j), and p t = pb + p m the total den¬ 
sity of monomers with pb = ^( r — R«j) be¬ 

ing the local density of brush monomers. (Here and in 
the following, the subscript b will be used to denote the 
brush chains, and to will denote the minority chain.) The 
first two terms in the Hamiltonian describe the Gaus¬ 
sian stretching energy (bonded interactions) between two 
neighboring monomers in the same chain. The third term 
represents the non-bonded effective interactions between 
polymer beads. For a good solvent the excluded volume 
parameter Ubare is larger than zero. The fourth term 
describes the adsorption between the substrate and the 
minority beads with a step-like adsorption potential, 


Hads(^) 


—e, 0 ^ z < a 
0, otherwise 


( 2 ) 


In the MC simulations, local densities are extracted 
from the position of the beads by using the Particle- 
to-Mesh technique I38|, which provides a way for the 
smoothing of density operators. In the present work we 
use the zeroth order scheme. The system is divided into 
cubic cells of size b 3 whose centers define the grid points, 
and the local densities evaluated by counting the total 
number of beads in the corresponding cells divided by 
the cell volume, are defined at these grid points. The 
smoothed monomer densities for the minority chain and 
brush chains are denoted as p m and pb, respectively. This 
implies that two beads interact only when they are in 
the same cell, and thus the size of the cell is indirectly 
linked to the interaction strength of non-bonded interac¬ 
tions. Higher order schemes are conceivable, but compu¬ 
tationally more expensive. After the density operators 
are smoothed over the cell volume, the Hamiltonian ac¬ 
quires the familiar form used in the SCF theory 31 fl . 
hence SCF calculations and MC simulations can be com¬ 
pared directly. 

In the present study, the brush chain length is fixed 
at Nb = 100. In most MC of calculations we choose 
the excluded volume parameter i>bare = 1 which would 
correspond to a Flory-Huggins parameter y = 0 in an 
explicit solvent solution (athermal solvent condition), 
and the size of the averaging cell is taken b = 1. We 
model relatively dense monodisperse polymer brushes 


in good solvent with the surface grafting density a = 
0.1, 0.2, 0.3, with the corresponding overlap parameters 
<jR 2 gb — 2, 3, 5, where R g b is the mean radius of gyration 
for an ideal brush polymer. The grafting points of the 
chains on the substrate were fixed on a regular square 
lattice.The system size was chosen L x = 40, L y = 40, 
and L z = 100 in most cases. To assess the influence of 
finite-size effects, we have also carried out simulations of 
larger systems with size L x = 60, L y = 60, and L z = 100 
for selected parameter values. The results were identical 
within the error. 

MC simulations [13 were carried out according to the 
standard Metropolis criterion. At every MC update, we 
try to move the position of one chosen monomer to a new 
position with a distance in space comparable to the bond 
length. This trial move results in an energy change in¬ 
cluding the bonded energy and non-bonded energy, and it 
is accepted or rejected according to the Boltzmann prob¬ 
ability in the usual way. In all the simulations 3 x 10 5 
MC steps per monomer were performed to equilibrate the 
system, and another 3 x 10 5 MC steps to extract statis¬ 
tical averages. In order to get good statistics sampling 
for 

the single minority chain, an additional 10 8 MC steps 
per monomer updating the minority chain monomers 
were performed during the process of evaluating statis¬ 
tical averages. Final statistical quantities were obtained 
by averaging the results from 48 separate independent 
MC runs. 


III. HOMOGENOUS BRUSH AND 
RENORMALIZATION OF THE SECOND VIRIAL 
COEFFICIENT 

We first analyze the simulation results for homogenous 
monodisperse brushes which are very well understood. In 
contrast to common MC polymer models our model con¬ 
tains an additional free parameter, the cell size b , defining 
the averaging volume for the local density operator, and 
the excluded volume parameter Ubare- Fig. [L] shows the 
density profiles of homogeneous brushes with two differ¬ 
ent grafting densities a = 0.1 and cr = 0.3, and various 
values of the averaging cell size b. It is clear that the 
choice of cell size has a pronounced effect on the brush 
density profile. To rationalize this effect we turn to the 
analytical SCF theory developed over 20 years ago 0 - 
0| and subsequently verified by simulations and experi¬ 
ments. 

Within the second virial approximation for good sol¬ 
vent conditions and assuming Gaussian elasticity for the 
chains, the density profile Pb(z) has a parabolic form 

«W-5S? (ff2 - z2> <3) 

where the brush height is H = (4n e fftr/7r 2 ) 1 / 3 A r f ) , and 
an effective interaction parameter. Gaussian elasticity is 
an important ingredient in the analytical theory and the 
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FIG. 1: Monomer density profiles for pure brushes at 
grafting density a = 0.1(a) and a = 0.3(b) from MC 
simulations obtained with interaction cells of different 
size b as indicated. 


fact that the MC simulation model is also based on a 
Hamiltonian with a Gaussian term describing chain con¬ 
nectivity allows a more direct comparison. It is known 
that corrections due to finite chain extensibility become 
important for dense grafting, but this is outside the scope 
of the present paper. On the other hand, modifications of 
the chain elasticity due to the excluded volume effects are 
naturally taken into account by our MC simulations. The 
mean-field potential associated with the density, profile, 
Vb{z) = v e spb(z) is also parabolic 


3tt 2 z 2 

m = v °-^Nf ( 4 ) 

where Vo = §(7ru e ffcr/2) 2//3 is the potential at the graft¬ 
ing surface. According to the theory (see for example 
[45ip. the brush density profiles evaluated for the same 
model at different values of a and Nf, should collapse 
in the rescaled coordinates vs a i/ z 3Nb with one ad¬ 
justable parameter related to the interaction parameter 
K = H/(a 1 / 3 N b ) = (4u e ff/7r 2 ) 1/3 


p{z) = JL 

o’ 2 / 3 2 K 3 


rV3jv h 


2 ' 


( 5 ) 


FIG. 2: Density profile of polymer brushes (a) in 
rescaled coordinates pb{z)/a 2 / 3 vs z/(a l ^ 3 Nb) 
calculated by MC simulations with Ubare = 1 for 
Nb = 50,100, 200 and different grafting densities 
a = 0.1, 0.2, 0.3. Dashed line shows the best fit to Eq. 

©, corresponding to K = 0.68. Panel (b) shows the 
dimensionless ratio S e ff/-Bbare as a function of b/a with 
B = i& 3 ( 1 — e~ v / b ) for different Ubare- These points are 
calculated from K obtained from the best fitting curves 
at fixed chain length Nb = 100. 


The rescaled MC density profiles with a = 0.1, 0.2, 0.3 
and Nb = 50, 100, 200 are presented in Fig. |2(a)| for 
a specific choice of the averaging cell size, 6 = 1 and 
the “bare” excluded volume parameter Ubare = 1 • The 
quality of the collapse is good, and one can extract the 
value of the effective interaction parameter v e ff = 0.776, 
which differs from Ubare- The same procedure was re¬ 
peated for several other values of the b parameter, and 
different bare excluded volume parameters fbare ■ To eval¬ 
uate the difference between v e g and Ubare, we calculate 
the corresponding exact second exact virial coefficient in 
the simulation model, taking into account the finite grid 
size b. For given interaction parameter v, it is given by 
B = \b 3 {\ — e~ v / b3 ). The ratio B e s/Bi, aie for the best fit 
values of are plotted as a function of ^ in Fig. |2(b)| 
The data for different values of Ubare collapse nicely. 

We conclude that the renormalization effect can be pre- 
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sented in terms of a dimensionless ratio £ eff , which is a 

-£>bare 

universal function of a single dimensionless parameter ^ 
(or, equivalently, by/N / R g ). Universality in this context 
does not mean a broad class of different models but rather 
that other model parameters, such as a and Nf, do not 
enter explicitly. It is clear from Fig. |2(b)~| that at cell size 
b/a = 2 or larger, the curve saturates and v e s approaches 
the bare value. Hence the system approaches mean-field 
behavior if the interaction range b of the monomers be¬ 
comes large, as one would expect. At smaller b/a , local 
monomer correlations become important and renormalize 
the effective interaction parameter. Renormalization of 
monomer interaction parameters has also been observed 
in Edwards-type simulations of polymer melts [38}. 


IV. MINORITY CHAIN IN A BRUSH: 
THEORETICAL BACKGROUND 


The theory was briefly sketched in Ref. E3- In this sec¬ 
tion, we present it in more detail. We consider a relatively 
dense monodisperse polymer brush in a good solvent con¬ 
taining a single surface-active minority chain. The length 
of the minority chain is denoted as N, while the length 
of the brush chain is denoted as Nf,. The minority chain 
and the brush chains are chemically different and have 
different interactions with the substrate, i.e., the sub¬ 
strate adsorbs minority chain monomers with strength e, 
while it is neutral to the brush chain monomers. Both 
the minority and the majority chains are flexible with 
the same statistical Kuhn length a which is taken as the 
unit length; the excluded volume parameter is also the 
same for both chain types and is positive corresponding 
to good solvent conditions. This is fully consistent with 
the MC model. 

We treat the inter- and intrachain interactions in the 
mean-field approximation and neglect the back effect of 
the change in the minority chain conformation on the sur¬ 
rounding brush. Hence, the conformation of a minority 
chain is affected by a fixed mean-field potential profile 
consisting of the repulsive contribution determined by 
the brush density and a short-ranged attraction due to 
the solid substrate. The minority chain itself is described 
by an ideal continuum model, since intrachain excluded 
volume effects are screened out considerably within the 
brush thickness. 

The analytical description of the minority polymer 
chain is based on the continuum approach. The Green’s 
function G(z,N), i.e. the total statistical weight of the 
minority chain grafted at the substrate (zg = 0) as a 
function of the free end position, z, is a solution of the 
Edwards equation 


dG(z, s) 
ds 


1 d 2 G{z,s) 
6 dz 2 


V(z)G(z,s) 


( 6 ) 


is the adsorption attraction potential, and 14 is the mean- 
field brush potential given by Eq. ©■ The adsorption 
potential I4ds(z) is described in the analytical SCF ap¬ 
proach as an attractive pseudopotential T4d 8 (z) oc —cS(z) 
where c is the adsorption interaction parameter in the 
continual model. This parameter was introduced by 
de Gennes [46} to replace the real adsorption potential 

through the boundary condition [G(z, s)] -1 dG g Z z ’ s ^ |_ Q = 

—c. Close to the substrate, z -C H, the brush po¬ 
tential changes very little, Vb(z) ~ Vo = kct 2 / 3 where 
K = I (■ nv e g/2 ) 2//3 ~ 1.71 is a shorthand notation for the 
numerical coefficient 

The adsorbed state is thus approximately described by 
the standard Green’s function modified by the potential 

V 0 0 

G ads (z) = 2 ce M-U,+c 2 /6)-«_ ( 7 ) 


Integration over z gives the partition function 

Q ads = 2e Ar( - y ° +c2 / 6) (8) 

In order to connect this expression to the MC model we 
note that c 2 /6 has the meaning of the negative free en¬ 
ergy of adsorption per segment (the chemical potential) 
in an asymptotically long chain, i.e., fi = c 2 /6. When 
comparing the actual brush density profile with the the¬ 
oretical parabolic shape, one concludes that the pseu¬ 
dopotential must account for the combined effect of the 
actual adsorption potential (a step of unit width with en¬ 
ergy e < 0) and the short-range depression in the brush 
density near the wall. It is known that for weak adsorp¬ 
tion, the free energy per monomer is determined by the 
crossover exponent </> : fi ~ (e—Sc) 1 ^. For an ideal chain, 
(j> = 1/2 while the value for a chain with excluded volume 
was a subject of extensive investigations and prolonged 
debates [H, |4t|. For any practical purposes, the adsorp¬ 
tion of relatively short chains is very accurately described 
by the ideal value <f> = 1/2. Here, for weak adsorption, 
we adopt the expression 

/r(e) = a(e - e c ) 2 (9) 


where a and e c are model-dependent constants to be 
extracted from the MC data. It is also known that at 
stronger adsorption, the free energy deviates from Eq. 
©, and for very large e, it can simply be estimated as 
fi = £ — const, where the constant shift is given by the 
limiting entropy difference per monomer between the coil 
and the fully adsorbed state. 

The state with the free end exposed at the brush edge 
has no contacts with the surface. Hence the Green’s func¬ 
tion is given by the known solution of the Edwards equa¬ 
tion for a purely parabolic potential and a neutral solid 
surface [la ] 


G ex (z,N) 


7r 
2 


3/2 


3 


Nb sin 



ze 


Z 


2 


taken at s = N, with the initial condition G(z, 0) = S(z). 
The total potential V(z) = Vb{z) + V ads {z), where U a d s (z) 


3tt f ( tvN A 
4 N b COt V 2 N b ) 


( 10 ) 
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For longer minority chains, N > Nb , the Green’s func¬ 
tion of Eq. ca shows unlimited monotonic increase with 
z since the solution refers to a potential extending to in¬ 
finity and to an infinitely extendable Gaussin chain. In a 
realistic situation the minority chain does not make ex¬ 
cursions well beyond z = H. To account for this fact we 
truncate G ex (z, N ) at z = H. If minority chains are just 
slightly longer A = N — Nb -C Nb, (which is the situation 
of interest in this paper), the Green’s function simplifies 
to 

G ex (z,N) = \ f(&) (n) 

{ 0, z>H 

This gives the partition function 

Qex = ^( 3N b) 1/2eV ° A ( 12 ) 

assuming e v ° A 1. Eqs (J8J) and o. clearly show that 
in the limit of Nb —> oo, a = const, N ^ b = const the 
free energies of the two states are extensive and switching 
become a classical first-order phase transition. 

Based on this analytical approach, the transition prop¬ 
erties between the adsorbed state and exposed state can 
be investigated. For example, one could define the tran¬ 
sition point as the condition Q a ds(c*) = Q e x- Omitting 
logarithmic correlations this gives fj* =Voj$. Under the 
condition (7^ C 1 one can use the asymptotic expression 
Eq. & t° obtain the adsorption energy at the transition 

-e* =-£ c +(/t/a) 1/2 cr 1/3 (-^) 7 (13) 

For larger grafting densities and chain length differences 
one expects deviation from the scaling form of of Eq. 
(fldl) . since the quadratic dependence of /z(e) may have 
a limited range. In the extreme case of a very long mi¬ 
nority chain, N Nb, the transition point satisfies the 
condition fi* = Vq- Since the mean-field potential Vq is 
in the range Vo ~ 0.36 — 0.8, the limiting value of e* 
remains of order 1. 

The sharpness of the transition is characterized by the 
width of the adsorption parameter change, Se, sufficient 
to produce a reliable switching from the adsorbed state 
with the average distance of the free end of the minority 
chain Z en <C H to the exposed state with Z en( j ~ H: 
Se = (jj dZ ^ d ) -1 where the derivative is evaluated at 
the transition point itself. A standard recipe for a two- 
state model [5Cj gives Se = (V dln Q^ )-i | _ i ) _1 which 
combined with Eq. ((5J results in 

Se = — {n/a) 1 ^ 2 a~ 1 ^{N/N) -1 ^ 2 (14) 

K 

Finally, the switching time is essentially determined 
by the barrier separating the two relevant states. Hence 
we estimate the barrier height at the transition point, 
Ubarrier, in terms of the change in the Green’s function of 
the exposed state, [/barrier = In ,, . it follows 

'- T ex(,^ — Z barrier y 


from Eq. HID that up to a numerical coefficient close to 
unity 

[/barrier ~ U 0 A ~ K( T 2 / 3 A (15) 


V. SIMULATION RESULTS AND DISCUSSION 

A. Adsorbed and exposed conformations of the 
minority chain 

Fig. [3] shows typical simulation snapshots of a brush 
containing a minority chain N = 120 (monomers) in an 
exposed and an adsorbed state. There is a substrate sur¬ 
face below and a virtual surface above showing the brush 
height H. As can be seen a small fraction of monomer 
units are located above H. These monomers form a tenu¬ 
ous brush exterior. The minority chain is shown in green 
color. It is longer than the majority brush chains with 
Nb = 100. We characterize the length of the minority 
chain by a parameter A = N — Nb (In the snapshot, A 
is given by A = 20). The minority chain in the exposed 
state has a stem-crown-like conformation. 

The main features of the snapshot can be also seen 
in Fig. [|] where the longitudinal monomer density pro¬ 
files are shown for minority chains with different A in 
the exposed state at e = 0 (panel a), and in the ad¬ 
sorbed state at £ = — 0.7(panel b). The density profile of 
a homogeneous brush is also shown, and for comparison 
both profiles are normalized in the same way. It is clear 
that at e = 0 and A = 0 the minority chain is identi¬ 
cal to brush chains. For longer minority chains, A > 0, 
the density profile deviates from a parabolic, monoton- 
ically decreasing form (apart from a narrow depletion 
layer near the substrate) shape. It is well understood 
that the standard parabolic shape is due to a specific 
very broad end-monomer distribution. “Partial” density 
profiles produced by chains with a given end position are 
initially increasing with z with a subsequent sharp drop. 
Only after averaging over end positions one obtains the 
actual density profile. The most pronounced effect is 
observed for the longest minority chain with A = 20. 
The profile can be understood if we take chain conforma¬ 
tions with the end-monomer positioned on average close 
to z = H (shown by the vertical dashed line) with a 
typical fluctuation of a few monomer lengths (~ 5). For 
shorter minority chains the average end-hight is below H 
and the fluctuations are relatively larger, which makes 
the difference to the majority chains less dramatic. In 
the strongly adsorbed state the density profile exponen¬ 
tially decreases with z and is insensitive to A, see Fig. 

M 

Even though the minority chains are highly dilute, they 
perturb the brush locally. Furthermore, the minority 
chain itself of course has lateral structure. Fig. [S] shows 
the lateral distribution of the end monomer P xy (r) as a 
function of the in-plane radial distance r = \jx 1 + y 2 
from the grafting point for different adsorption energies 
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(b) 


FIG. 3: Simulation snapshots of a polymer brush with 
Nt = 100 and grafting density a = 0.1 containing a 
minority chain with N = 120 in the exposed state (a) at 
e = 0 and the adsorbed state (b) at e = —0.7. The 
minority chain beads are colored green, while the brush 
chain beads are colored red and made transparent. 
Brush height H is shown by empty square above for 
orientation. 


e at grafting density a = 0.2 and chain length N = 120. 
Here P xy (r) was obtained by counting how often the mi¬ 
nority chain end was located in a cylindrical shell rang¬ 
ing from r to r + A?’ and dividing this by the volume 
of the shell and the total sampling number. The thick¬ 
ness of the shell was chosen as A r = 1. The results 
for £ = 0 in Fig. [5] show that P xy (r ) exhibits Gaussian 
behavior for nonadsorbing chains. With increasing ad¬ 
sorption strength, the distribution becomes broader, and 
a depletion region develops close to r = 0: Chain ends are 
pushed away from the area close to the grafting point, due 
to the fact that this area is already covered by adsorbed 
monomers. This is also shown by the radial distribution 
density P r (r) = 27rrP a ,j / (r), which gives the probability 
to find the free end of the minority chain located at dis¬ 
tances ranging from r to r + A r with A r = 1. At strong 
enough adsorption the maximum of P r (r) moves towards 
larger values of r since the flattened adsorbed state dom- 



z 



z 

(b) 

FIG. 4: Longitudinal monomer density profiles of 
minority chains of different length N = Ny + A in the 
exposed state at e = 0.0 (a), and in the adsorbed state 
at £ = —0.7 (b) as obtained from MC simulations. The 
brush parameters are Ni, = 100 at cr = 0.1. The density 
profile of the brush chain normalised to unity with is 
shown for comparison by solid red line. Brush height 
H « 32 is shown by vertical dashed line. 


mates (Fig. O inset). 

To measure more quantitatively how the minority 
chain extends in the lateral direction, we calculate the 
parallel root-mean-squared radius of gyration defined as 

i?H = yjRg X + Rg y , where R gx and Rgy are the mean- 
square radius of gyration of the minority chain in the x 
and y directions, separately. 

Fig.[5]shows i?.|| as a function of the adsorption strength 
at cr = 0.1 for different minority chain lengths. From 
these curves it can be seen that the lateral size of the mi¬ 
nority chain remains almost unchanged up to £ ~ —0.4, 
and then starts to increase, which is another signature of 
the conformational transition from the desorbed to the 
adsorbed state. Fig.[7]shows the lateral size of the minor¬ 
ity chain in the exposed state with £ = 0 and in a strongly 
adsorbed state with —£ = 0.7, as a function of the mi¬ 
nority chain length for several grafting densities. In both 
the exposed state and the adsorbed state, P|| ( N ) is con¬ 
sistent with a power law. The best fit for the exposed 



























r 

FIG. 5: Lateral distribution of the free end of the 
minority chain as a function of radial (in-plane) 
distance from the grafting point r for a = 0.2, N = 120, 
and different adsorption strengths as indicated. The 
inset shows the lateral probability given by 
P r = 2TrrP xy (r)Ar , where Ar=l. The data are 
obtained by MC simulations 



FIG. 6: Lateral size of the minority chain defined as 
.R|| = yjFt-lx + Rgy vs the adsorption strength for 
a = 0.10 and different minority chain lengths. The data 
are obtained by MC simulations 

state (e = 0) gives R\\ oc N Ve , with v e « 0.55, which is 
larger than the Flory exponent expected for ideal chains 
(v = 0.5), and smaller than that expected for free chains 
in a good solvent {v = 0.588). This can be understood 
since part of the chain well within the brush experiences 
semi-dilute conditions with screened excluded volume in¬ 
teractions while the part near the tenuous brush exterior 
is swollen. As one increases the grafting density at fixed 
length, the lateral size of the exposed chain decreases. 
This has two reasons: First, the brush thickens with in¬ 
creasing grafting density, and consequently the fraction of 
swollen minority chain part outside the brush (the crown) 
decreases. Second, the monomer density inside the (semi¬ 
dilute) brush increases, which leads to a stronger screen¬ 
ing of excluded volume interactions and a decrease of 


lateral swelling of the inner part of the minority chain 
(the stem). According to Fig. [Gl the shrinking does not 
depend noticeably on the chain length A = N — A), be¬ 
tween minority chain and brush chains, i.e. , on the size 
of the crown. Hence we conclude that the second effect 
is probably dominant. 

For strongly adsorbed chains (—e = 0.7), we find a 
relation R\\ oc N Va , with v a w 0.73, which is close to the 
Flory exponent for two dimensional self-avoiding chains 
{y = 0.75). In the strongly adsorbed state the lateral 
extension is almost insensitive to the grafting density, 
see Fig. [7J upper line. From the scaling point of view 
the adsorption blob size is quite small and thus screening 
of the excluded volume interactions by other chains is 
ineffective. 



FIG. 7: Lateral size of the minority chain vs chain 
length for different grafting density in the exposed state 
(e = 0) and in a strongly adsorbed state —e = 0.7. The 
symbols denote the results from MC simulations, lines 
are obtained by a linear fitting after taking the 
logarithm of the original data 


B. Transition properties 

Fig. |8(a)| presents the average distance Z en d from the 
free end of the adsorption-active minority chains to the 
solid substrate as a function of the adsorption strength 
—e, for minority chains of several lengths in a brush with 
Nt = 100 at a = 0.2. The two lower curves represent 
the behavior of shorter minority chains with N = 50 
and N = 90. In the desorbed state, short minority 
chains exist in a coil-like conformation and hence the ad¬ 
sorption transition is expected to be continuous. The 
curves are further smeared out due to the finite length of 
the chain. In contrast, even a small positive increment 
A = N — Nb in the length of minority chain results in a 
well-pronounced sharp transition from an exposed state 
in which the free end of the minority chain is localized at 
the outer space of the brush to an adsorbed state where 
the chain is localized close to the substrate as illustrated 
in Fig. |T] 
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Fig. |8(b)~| compares the MC results with the prediction 
of the three-dimensional SCF, using the interaction pa¬ 
rameter u e fi extracted from the MC simulations of pure 
brushes (Fig. |2(bJ| ) . The SCF results differ from the MC 
date in two important aspects. First, they do no capture 
the attractive contact interaction of the minority chain 
with the surface accurately. This is not surprising, given 
that the interaction vanishes on the scale of one grid size, 
and it leads to a shift of the adsorption curves by a con¬ 
stant value £o ~ 0.1. Second, the density fluctuations 
at the brush surface (see Fig. |2(a)[ ) effectively reduce the 
(spatial) range of the repulsive potential, hence the mi¬ 
nority chain is less stretched in the exposed state. This 
effect is even observed for minority chains that have the 
same length than the brush chains, N = Nb. We will 
see below that this leads to a reduction of the free en¬ 
ergy barrier between the exposed and the adsorbed state. 
After shifting the SCF curves by £q and rescaling them 
with a factor a(N) which was obtained empirically, they 
are in very good agreement with the MC results. 

Curves such as those shown in Fig. [8j which give the 
average distance of the minority chain free end from the 
surface as a function of the adsorption strength, can be 
used to localize the crossover between the adsorbed and 
the desorbed state. The transition point — e* is then de¬ 
fined as the point where the slope of the curve is maximal 
(for A > 0). The value of the slope at this point is used 
to extract the width of the transition, 5e. This is done as 
follows. First we find the point where Z eTL d vs —£ has the 
maximum slope, and calculate this slope as l. Then we 
draw a line through this maximum slope point with the 
slope l and find its intersection with (i) the abscissa and 
(ii) with a line parallel to the abscissa that corresponds 
to the value of Z en d at £ = 0. The absolute value of 
the difference of the adsorption strength at these two in¬ 
tersection points (i.e., difference of their abscissa values) 
defines the transition width 


8e = 


•^end|e=0 


l 


(16) 


In order to study the transition mechanism we analyze 
the distribution of minority chain ends, which we denote 
as P z . The product P z dz represents the probability to 
find the free end in a layer with its z coordinate ranging 
from z to z + dz. In the MC scheme, P z is evaluated as 
follows: We split the simulation box along the z direction 
into n z layers, where each layer has the volume L x * L y * 
L z /n z and an index k ranging from 0 to n z — 1. Then 
P z {k) is obtained by counting how often the free end was 
located in the fcth layer, divided by dz and by the total 
size of the sample. 

Fig. 9(a) shows examples of distributions P z (z) for a 
minority chain of length N = 110 in a brush with grafting 
density a = 0.3. One can clearly see that two popula¬ 
tions of different conformations (adsorbed and exposed) 
are involved in the transition. In the transition region 
these populations coexist which is reflected by a bimodal 
distribution. Specifically one can see that at low adsorp¬ 
tion strength (e.g., £ = —0.26), the chain end is mostly 
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FIG. 8: Average distance Z en( j of the free end of the 
minority chain from the surface vs adsorption strength 
for different chain lengths with grafting density a = 0.2 
obtained from MC simulations (a). Superposition of a 
subset of these curves (MC) with the corresponding 
curves obtained by 3d SCF calculations (SCF), where 
the SCF curves have been shifted by £ 0 = 0.1 and 
rescaled by an empirical factor a(N), given in the inset. 
The SCF calculations were done with the renormalized 
interaction parameter u e ff extracted from Fig. |2(b)[ 


located outside of the brush, whereas at high adsorp¬ 
tion strength (e.g., £ = —0.6), the chain end is local¬ 
ized close to the surface. At intermediate adsorption 
strengths (e ~ —0.42), the adsorbed and the exposed 
state coexist and the minority chain switches back and 
forth from one to the other. 

A conventional way to analyze a phase transition is to 
construct a “Landau free energy” from the logarithm of 
the statistical distribution of the order parameter a 
We will use the position z of the free end of the minority 
chain as the order parameter. In this case the Landau 
free energy is — In P z {z). Selected free energy curves are 
shown in Fig. |9(b)] The free energy as a function of chain 
end position has one minimum close to 2 = 0 for large 
£ (adsorbed state), and one minimum close to 2 = 40 
for small £ (exposed state). For intermediate £, these 
two minima compete with each other. An alternative 
definition of the transition point e* is the point where 
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(a) 



(b) 


FIG. 9: Longitudinal distribution of the free end of the 
minority chain obtained from MC simulations as 
functions of z at cr = 0.3, A = 10, for different 
adsorption strength (a). The Landau free energy is 
obtained by taking the negative logarithm of the 
distribution P z in Fig. |9(a)] The black dashed vertical 
line is shown to denote the brush height H m 46 at 
N b = 100, cr = 0.3. 


both minima have equal depth. However, we wish to 
stress that at finite chain length, there is no classical 
phase transition, but rather a crossover point between 
two regimes. 

The height of the free energy barrier between the ex¬ 
posed state and the adsorbed state, C/barrien is extracted 
from the Landau function curve at this transition point. 
The equal depth condition for the transition point makes 
the barrier height definition unambiguous. In the follow¬ 
ing we compare the transition characteristics as obtained 
by simulations to the scaling prediction of the analytical 
theory laid out in Section [TV] 

Fig. [TO] displays the transition point, i.e. the adsorp¬ 
tion energy —e* as a function of the scaling parameter 

cr 1 / 3 (t|) combining the grafting density, and the ratio 
of the minority-to-majority chain lengths. Two defini¬ 
tions of —e* have been proposed above, and data points 
corresponding to both definitions are presented. Each set 
collapses into a single master curve close to the straight 
line suggested by Eq. m but due to the effect of finite 



FIG. 10: Transition points as a function of the scaling 
parameter o 1 ^ 2 (A/N) 1 ^ 2 obtained from finding the 
value e* corresponding to the criterion of equal height 
in P z (see Fig. |9(a)| ) (filled symbols), o r to t he criterion 
of maximum slope of Z en( j (see Fig. 8(a) I (empty 
symbols). The solid and dash lines are used to guide 
the eyes and they have almost the same slope 


chain length the lines are separate. We expect the differ¬ 
ence to vanish in the appropriate thermodynamic limit 
N —> oo, N b —> oo, ^ = const, o = const, where the tran¬ 
sition becomes truly sharp. The general trend is that the 
transition point shifts to larger adsorption strength |e*| 
with increasing minority chain length and grafting den¬ 
sity. This is to be expected, since both factors stabilize 
the exposed state against the competing adsorbed state. 

In Fig. [TOl the intercept and the slope are linked to the 
model-specific dependence of the adsorbed chain’s free 
energy on e, which were therefore treated as free param¬ 
eters by the theory. One has to be cautious, however, 
not to overestimate the agreement between theory and 
simulations. The scaling predictions rely on two assump¬ 
tions: (i) logarithmic corrections are neglected, implying 
that the relevant free energies per chain are much larger 
than 1 (in kT units), Vo A 1; (ii) the adsorption is 
assumed to be weak, such that the expansion of fiie) 
in Eq. fliJ]) can be applied, -C 1. At brush chain 
lengths N « N b = 100, it is difficult to satisfy both con¬ 
ditions simultaneously. In systems with small values of 
A and/or small o, the first condition becomes question¬ 
able, whereas in systems with larger A and/or larger a, 
the second condition may be violated. Therefore, the 
observed reasonably good collapse onto a straight line 
might be the result of cancellation effects. 

The two most important characteristics of the tran¬ 
sition, as far as potential applications are concerned, 
are the transition sharpness (or effective width) and the 
waiting time characterizing the transition kinetics. The 
present work does not include dynamic simulations, so 
we use the free energy barrier at the transition as an 
indirect measure of the waiting time. The transition 
width, denoted as 6e and charac terized by the inverse 
slope of Z e nd(£), is plotted in Fig. 11(a) as a function of 
the minority chain length increment A, while the barrier 
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FIG. 11: Transition sharpness (a) and barrier height 
between coexisting states (b) as a function of the chain 
length difference A = N — Nf, between minority chain 
and brush chains for different grafting densities as 
indicated, according to SCF theory (lines) and MC 
simulations (symbols). 


height (in units of kT) is displayed in the Fig. |ll(b) We 
present both the MC results (symbols) and the SCF re¬ 
sults (lines). The transition becomes sharper for longer 
minority chains and larger grafting densities. The differ¬ 
ence between the MC and the SCF results for Se (Fig. 
11 (a)|) i s noticeable but not dramatic. As discussed in 
Sec.lIVIfEa. itTUl b the transition sharpness is essentially 
determined by the response of the free energy of the ad¬ 
sorbed state to variations of the adsorption strength e, 
which does not seem to be affected by fluctuations very 
much. In contrast, the effect of fluctuations on the free 
energy barrier is quite strong. The free energy barrier 
grows approximately linearly with A and also increases 
with the grafting density. The absolute values of the bar¬ 
rier height obtained by MC simulations are smaller than 
those from SCF method by roughly a factor of 2 (Fig. 
|ll(b)[ ) . This is clearly a pronounced effect of density fluc¬ 
tuations implying strong consequences on the transition 
kinetics. 

The same data are replotted in Fig. m with scal¬ 
ing coordinates as suggested by the analytical theory, see 
Eqs (TUO) . (1151) . The curve for Se vs. er -1 / 3 (IVA) -1 ^ 2 is 


FIG. 12: Same data as Fig. [TTJ shown in the scaling 
plot suggested by the analytical predictions in Sec. HVl 
with the scaling variable (cr 2//3 A). 


expected to pass through the origin since in the thermo¬ 
dynamic limit, one has cr -1 / 3 (JVA) -1 ^ 2 —> 0, and the 
transition becomes jump-like with Se —> 0. At finite 
brush chain length Nf, and for larger values of the scal¬ 
ing parameter er -1 / 3 (_/VA) -1 ^ 2 finite chain length effects 
and logarithmic corrections become important. This pre¬ 
sumably explains why the data do not collapse on a 
master curve, and why the SCF data don’t even seem 
to converge towards the expected behavior Se —> 0 at 
er- 1 / 3 (fVA)” 1/2 -s. 0. The MC data are better com¬ 
patible with a limiting straight line passing through the 
origin, and the slope has the correct order of magnitude. 
According to the theoretical expectation (Eqs (fT3l) and 
(HI) the slope in Fig. |12(a)| should exceed that in Fig. 
EH by a factor of ^ ~ 1.2. The actual slopes differ by a 
factor of 2.3. 

Fig. |12(bJ| displays the energy barriers at the transition 
as a function of cr 2//3 A. The SCF data are quite close 
to collapsing onto a single straight line, as predicted by 
Eq. (1151) while the MC data are not. Although the di¬ 
rect proportionality to the chain length increment A is 
still satisfied reasonably well, the effect of increasing the 
grafting density is stronger than predicted by the theory. 
This is probably due to the fact that the values of the 
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barrier height f/barrier obtained from MC simulations he 
in the range between 0.5 and 5 whereas the theory 
assumes that the two states are well separated, imply¬ 
ing C/barrier/fcsr 3> 1. Thus the asymptotic behavior 
assumed by the theory is not yet reached, which may 
explain the poor data collapse. 

It is worth noting that the data do seem to obey an 
apparent scaling f/barrier ~ crA as illustrated by Fig. Q1J] 
We are not aware of a theoretical reason for such a scal¬ 
ing. The observed apparent scaling might be the purely 
fortuitous result of various corrections to the analytical 
theory, which relies on several somewhat questionable as¬ 
sumptions as discussed above. 




FIG. 13: Same data of the transition barrier as Fig. [T21 
scaled as suggested by a scaling hypothesis with 
modified scaling variable a A. 


VI. CONCLUSION AND REMARKS 

To summarize, we have investigated in detail con¬ 
formational properties and phase transition behavior of 
a type of stimuli-responsive polymer materials using 
particle-based Monte Carlo simulations and continuum- 
based 3-dimensional self-consistent field theory. This 
type of stimuli-responsive materials is designed based 
on a polymer brush with a small fraction of adsorption- 
active minority chains introduced in the brush, and the 
minority chain can work as a responsive sensor. The 
basic mechanism for this responsive sensor relies on the 
conformational switch of the minority chain. The adsorp¬ 
tion between the substrate and the minority chain serves 
as the trigger that switches the chain conformations. In 
practice, the trigger could also be a temperature or sol¬ 
vent composition change. Upon varying the adsorption 
strength, two states were observed, one being an exposed 
state, in which the free end of the minority chain is lo¬ 
cated at the outside surface of the brush and possesses 
a stem-crown-like configuration, and the other being an 
adsorbed state, in which the minority chain is located at 
the adsorbing substrate in a nearly 2-dimensional con¬ 
finement. 

One important purpose of the present study was to 


FIG. 14: Barrier height versus sharpness of the 
transition according to SCF calculations (lines) and MC 
simulations (symbols). At given sharpness Se, the 
barrier height is significantly smaller in the MC 
simulations. 


highlight the effect of fluctuations by comparing two nu¬ 
merical methods applied to one and the same model, 
namely MC simulations and 3d SCF calculations. The 
comparison shows that the SCF theory roughly captures 
the sharpness of the transition ( Fig. |12(a) but overes¬ 
timates the barrier height (Fig. 12(b)[ ). As discussed 
earlier, densit y fluctuations at the outer brush surface 
(see Fig. |2(a)[ ) effectively reduce the range of the repul¬ 
sive barrier created by the brush and hence the degree of 
chain stretching in the exposed state (see Fig. |8(b)j ). As 
a consequence, the free energy barrier between the ad¬ 
sorbed and the exposed state is also reduced. The prac¬ 
tical importance of this effect can be seen when plotting 
the barrier height against the sharpness of the transi¬ 
tion (Fig. [14]) . The transition sharpness is a measure for 
the sensitivity of the switch to small changes in the en¬ 
vironment. It can be tuned by adjusting the minority 
chain length or the grafting density. Assuming that it 
has been tuned to a certain value, Fig. [14] demonstrates 
that the corresponding barrier height between coexist¬ 
ing states at the transition is greatly reduced in the MC 
simulations compared to the SCF prediction. Since the 
barrier height determines the time required for switching 
between states, this can drastically reduce the response 
times of the switch. For example, for Se = 0.12, the re¬ 
duction can be as high as 7fcsT, and according to a sim¬ 
ple Arrhenius-type estimate this would lead to a speedup 
of the switching time by three orders of magnitude. 


The switching time is one major issue for the quality 
of sensors or other responsive materials. In our previous 
work m, we performed overdamped Brownian dynam¬ 
ics simulations of a single switch chain in a static brush 
potential to estimate the switching time. By tracking the 
position of the free end bead of the minority chain, we 
found that this bead could jump between two well sepa¬ 
rated positions (which correspond to the exposed and ad¬ 
sorbed state, respectively) during a time about 10 6 r mono , 
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where r mono is the characteristic relaxation time of a sin¬ 
gle monomer, which is typically on the order of 10 -9 sec¬ 
onds for flexible chains. We thus estimated the switch¬ 
ing time to be on the order of milliseconds. However, in 
that calculation, the brush was represented in a simplified 
manner by a parabolic profile and the coupling between 
the minority chain and the brush chains was neglected, 
which presumably leads to a large deviation from the true 
switch time. Molecular dynamic simulations that model 
a brush-minority chain system in a dynamic manner are 
clearly desirable. 

In the present study, we restrict ourselves to systems 
with good solvent conditions, and assume that two-body 
interactions dominate, thus we only consider excluded 
volume interactions. However, for dense brushes, higher 
order contributions will become important. Higher or¬ 
der terms in the free energy virial expansion provide cor¬ 
rections that lead to some change in the brush height 
and in the overall free energy of the competing states. 
Hence one would expect corrections to the properties of 
the switching transition, in particular the barrier height 
at transition. In the case of very dense grafting these 
corrections may be significant but this falls outside the 
scope of the present paper. 

Besides fluctuations, another factor that must be taken 
into account in real systems is polydispersity. In the 
present study, the brush polymers were taken to be 
monodisperse. Real polymer brushes, however, are al¬ 
ways polydisperse, and it has been demonstrated that 
polydispersity significantly alters both the equilibrium 
properties [52], |53| and dynamical properties [5#l5q ] of 
materials. Using a one-dimensional SCF theory, we 
have performed a preliminary study of a system with 
an adsorption-active chain embedded in a polydisperse 
brush with a continuous Schulz-Zimm length distribu¬ 
tion. We found that for A > 15, the switch sharpness 
is similar to that of the monodisperse brush obtained by 
1-dimensional SCF theory; however, the switch barrier 
is reduced. This suggests that the switching mechanism 
should be at least as effective in polydisperse brushes 
than in monodisperse brushes, the switching might even 
be faster. Further and more detailed studies are currently 
under way. 

An important question pertaining to the MC scheme 
used in our work is the choice of the size of the density av¬ 
eraging cell, 5, which has the meaning of an interaction 
range. In section 1III1 we have shown that the effective 
excluded volume parameter appears as the result of a 
renormalization of the “bare” parameter due to density 
fluctuations. For a given set of parameters A/ and era 2 , 
where a is the statistical segment length unit, the brush 
profile is uniquely defined by u e ff or, more precisely, the 
corresponding second virial coefficient in the MC model, 
B P ff, while the renormalization factor t f eff was shown to 
be a function of a single dimensionless variable -. The 
renormalization of monomer interaction parameters has 
been studied for single flexible swollen chains [57|, in the 
context of dense polymer systems [13, and in the field- 


theoretical context [58M6l[. We would like to emphasize 
that polymer brushes simulated by MC in the Laradji 
version [2fil ] represent a class of model systems where 
renormalization effects are quite strong and can be ex¬ 
tracted easily. The basic tenet is that large-scale prop¬ 
erties should depend only on the renormalized excluded 
volume parameter while its bare value may become rel¬ 
evant on small length scales comparable to a or b. One 
consequence in our system is that the brush profiles very 
close to the outer surface cannot be described by a sin¬ 
gle renormalized u e ff (see Fig. |2(a)[ , and that the surface 
interaction parameter e has to be renormalized indepen¬ 
dently from Ubare (as has been done in Fig. |8(b)| ). 

However, even after taking into account the renormal¬ 
ization of interaction parameters, we find that the prop¬ 
erties of chains in close vicinity and above the brush sur¬ 
face are still strongly influenced by density fluctuations. 
The density fluctuations at the brush surface significantly 
reduce the free energy barrier between the exposed and 
the adsorbed state in Fig. [Td] Thus, the minority chain 
serves as a delicate probe of both local and large-scale 
fluctuation effects. 
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Appendix A: SCF theory 

The SCF method used in the present work is formu¬ 
lated in the continuous 3d space, starting from the same 
Hamiltonian with that in the MC method mentioned in 
the text. In the SCF approach |62[ - |64| ]. the partition 
function Z = (C{r} ex p(~ /3%({R})) is rewritten as a 
functional integral, i.e., Z = f 'Duj'Dpt exp(—/9J 7 ), with 
the free energy functional expressed as 

/3B= l -j dvpl (r) — J dru(r)p t (r) 

- rib In Qb[w\ - In Q m [to + C/ a ds], (Al) 

where pt is the total density distribution, u is the corre¬ 
sponding conjugate auxiliary potential, and Qb , Q m are 
the single chain partition functions. In the calculation we 
set the excluded volume parameter v = u e ff (where u e ff is 
obtained from the MC simulation). Extremizing the free 
energy functional with respect to the fluctuating fields, 
i.e., the density pt and the potential ui, one obtains the 
SCF equations 

w = vp t 
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rNb 

Pt=Pb + Pm= / dsql(r,s)q b (r, N b - s) 

Jo 

nN 

+ / dsql l (T,s)q m (T,N - s) (A2) 

Jo 

where qj^ and q y (/x = 6, m) satisfy the modified diffusion 
equation 

d a 2 

—g M (r, s) = —V 1 2 3 4 5 6 7 8 9 10 11 12 g M (r, s) - w M g(r, s), (A3) 

J^( r ’ s ) = y y2 ^( r ’ s ) “ w /^( r > s )> ( A4 ) 

with = oj for yi = b, = w + {/ a d s for fj, = to, 
Dirichlet boundary conditions along the z direction, and 
periodic boundary conditions along the x and y direc¬ 
tions. The functions q M and q t (/r = 6, m) correspond to 
spatial integrals of the single-chain propagator between 
monomer s and the free or grafted end, respectively, over 
all possible end positions. Thus the initial condition for 
q :4 (free end) is g M (r, 0) = 1, and the initial condition 
for qj^ (grafted end) is chosen self-consistently such that 
the density of graft points is reproduced correctly, i.e., 
ql( r, 0) qb(r,Nb) = aS(z — zq) for the brush chains, and 
^(r, 0) q m (r,N) = S(r — r 0 ) for the minority chain, 
where ro = (0,0, zq) is positioned at the center of the 
substrate. The transition properties are extracted from 
the corresponding minority chain propagator. For ex¬ 
ample, one can calculate the distribution of free end via 
Pend{z) = f da; dyq^lr, N)/ f drq} n (r,N). From these 
distributions the transition point and transition barrier 


can be obtained in a similar way as that used in the MC 
simulations. 

The SCF equations are closed and iterative meth¬ 
ods are used to find their solutions. The auxiliary po¬ 
tential ui is used as the iterative variable, and a new 
potential for the next iteration is generated a/ ra+1 ) = 
lo (") + \(vp[ n ^ — u^), where A is a small number control¬ 
ling the stability of the iteration and the speed of converg¬ 
ing. In most of the cases we use A = 0.03. The iteration 
process stops if the iteration error is smaller than 10 -5 . 
There are several ways to solve the modified diffusion 
equation numerically, for example, real space finite differ¬ 
ence schemes [gl], pseudo-spectral schemes [H, [§7|, and 
so forth. Usually, pseudo-spectral schemes are reliable 
and efficient. However, in our system, we encountered 
numerical problems when chains are strongly stretched 
(e.g. at large grafting densities, or when chains are long). 
In such cases, the propagators sometimes assumed neg¬ 
ative values, and the SCF equations did not converge. 
Similar problems were reported in Refs. [33l J34|. There¬ 
fore, we adopt the Doublas-Brian scheme |68j], a real- 
space finite-difference method which did not suffer from 
this problem. The volume of the simulation system is 
chosen equal to that of the MC method. It is divided 
uniformly into n x * n y * n z = 40 * 40 * 200 grid cells, 
and continuous quantities are discretized on the vortices 
of these cells. The chain contour including the minority 
chain is uniformly discretized into 400 steps. About 300 
hundred iteration steps are necessary for the iteration to 
reach convergence. 
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